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I.  INTRODUCTION 

When  an  energetic  atomic  projectile  travels  through  an  amorphous 
target,  it  will  normally  be  deflected  at  seme  angle  from  the  original 
direction  with  some  energy  loss.    The  energy  loss  is  due  to  two  factors; 
nuclear  energy  loss  (associated  with  elastic  scattering  from  the  target 
particles)  and  electronic  energy  loss  (associated  with  inelastic  scattering- 
from  the  target  particles).    The  elastic  scattering  between  two  atoms  is 
dependent  on  the  following  factors:     the  mass  of  the  projectile  and  of  the 
target  atom,  the  energy  of  the  projectile ,  and  the  interaction  potential 
between  the  projectile  and  the  target  atom. 

The  deflection  angle  and  the  differential  cross  section  can  be 
calculated  by  classical  theory1'2  and  by  quantum  theory.2^    In  the  classical 
theory  the  orbit  or  path  of  motion  and  the  momentum  of  the  projectile  can  be 
specified  at  every  point.    In  quantum  theory,  however,  the  momentum  and 
position  of  the  projectile  can  only  be  specified  consistent  with  the 
uncertainty  principle.    The  phase-shift  method  is  exact,  but  it  can  only 
be  used  conveniently  when  the  energy  of  the  incoming  particle  is  small. 
When  the  kinetic  energy  of  the  incoming  particle  is  large  compared  to  the 
scattering  potential  the  Born  approximation  introduces  little  error.  The 
classical  approximation  is  valid  when  the  following  two  conditions  are 
satisfied:    l)  the  deBroglie  wavelength  of  the  projectile,  X,  is  negligible 
compared  to  the  size  of  the  scattering  center,  and  2)  the  collision  is  well 
defined  within  the  limitations  of  the  uncertainty  principle .S>6  According 
to  Bohr ,7  the  second  condition  is  satisfied  for  the  screened  Bohr  potential 
when  the  scattering  angle,  G,  is  larger  than 

9*  «  X/2tra 


where  a  is  the  size  of  the  scattering  center  (a  =  ^pfo**  ZLV\\^'  ^ 
an  illustration  of  this  limitation,  consider  a  scattering  of  100-kev  oxygen 
ions  (Zx  =  8)  with  carbon  atoms  (Z2  =  6).    The  above  formula  gives  0*  =  0.25 
degrees  and  the  classical  approximation  is  valid  for  all  angles  greater  than 
0.2$  degrees.    On  the  other  hand,  the  Born  approximation  is  valid  at  small 
scattering  angles  (when  the  phase  shifts  are  small).    The  two  methods 
complement  each  other.    There  exists,  however,  a  small  range  of  scattering 
angles  where  neither  procedure  is  valid. 

SverhartS  has  calculated  the  scattering  angle  as  a  function  of  the 
impact  parameter  with  the  screened  Bohr  potential  (which  can  be  written  in 
closed  form).    These  calculations  are  based  on  the  classical  theory.  There 
exist  other  forms  of  atomic  potentials  e.g.  Thomas -Fermi,9' 10  Firsov,11  and 
Thomas-Fermi-Dirac12  potentials.    Moliere13'1"  suggested  that  the  Thomas- 
Fermi  function  could  be  approximated  by  a  sum  of  three  exponentials.  More 
recently  there  has  been  increased  interest  in  forms  of  atomic  potentials 
which  can  be  expressed  in  closed  form.1!?    Csavinszky1^  reports  the  modi- 
fication of  the  Firsov  potential  which  is  expressed  as  the  sum  of  two 
exponential  terms.    Roberts1?  proposes  an  approximation  of  the  Thomas-Fermi 
function  which  consists  of  a  single  exponential  term.18    An  empirical  fitting 
of  the  Firsov  potential  over  a  large  range  of  intemuclear  distance  (but 
not  valid  for  small  r)  is  recently  reported  by  Wedepohl.19    The  potentials  - 
of  Csavinszky,  Roberts,  and  V.'edepohl  have  mors  desirable  properties  at 
larpe  r.    We  note  that  the  scattering  calculations  between  atoms  even  with 
the  classical  model  are  not  very  extensive,  and  the  sensitivity  of  the 
numerical  results  is  not  investigated. 

It  is  the  motivation  of  this  work  to  investigate  the  dependence  of 
the  scattering  angle  and  the  differential  cross  section  (with  a  classical 


description  of  elastic  scattering)  for  ions  in  the  kev  energy  region  on 
different  forms  of  potentials.    Section  II  contains  the  main  points  of  the 
classical  theory,  and  section  III  describes  the  procedure  used  in  the 
present  calculations.    Numerical  results  are  presented  in  section  IV, 
followed  by  a  discussion. 


II.    CLASSICAL  THEORY  OF  SCATTERING 

In  the  laboratory  system  the  scattering  is  considered  from  a  stationary 
target.    Momentum  is  transferred  from  the  projectile  to  the  target  atom, 
which  recoils.    The  following  quantities  enter  the  calculations  in  a  two- 
body  problem:    px,  r-^  p2»  and  ^2  (the  momentum  and  Fosition  vectors  of  the 
projectile  and  of  the  tar-get  atom,  respectively).    However,  the  center  of 
mass  of  the  two  particles  is  considered  at  rest  in  the  center  of  mass  system. 
The  two-body  problem  can  be  reduced  to  an  one-body  problem  provided  the 
interaction  depends  only  on  the  relative  co-ordinate  r  =         -  r^).  In 
the  center  of  mass  system  the  co-ordinates        and  r2  are  replaced  by  the 
co-ordinate  of  the  center  of  mass  R  =  (mi?i  +  m2r2)/M  where  M  =  mi  +  m2 
and  by  the  co-ordinate  r.    The  motion  of  the  center  of  mass,  which  is  the 
motion  of  a  free  particle  traveling  at  a  constant  velocity  Vc  =  m-jv/14 
where  v  is  the  initial  velocity  of  the  projectile  in  the  laboratory  system, 
can  be  separated  out  and  one  has  an  equivalent  one-body  problem  with  an 
effective  mass  m  =  m^/M.    The  experiments  are  done  in  the  laboratory 
system  end  therefore  we  note  below  the  relevant  relations  to  transform  the 
various  quantities  from  the  center  of  mass  system  to  the  laboratory  system. 
The  relation  between  the  scattering  angles  in  the  laboratory  system,  V  , 
and  in  the  center  of  mass  system,  0,  is 


tan  V  5 


sin  e  (i) 
cos  &  +  ~6 


X    is  defined  as  (BLXSx.  .  — £ — "/where  Q  is  the  amount  of  energy  converted 
from  internal  energy  to  kinetic  energy  of  the  emergent  particle.  For 
elastic  scattering,  Q  =  0,  m1  =  m^,  m2  =  m^,  ^  reduces  to  n^/ny  The 
initial  kinetic  energy  in  the  center  of  mass  system  is 


S  =  m2Eo/H  (2) 
where  EQ  is  the  initial  energy  of  the  projectile  in  the  laboratory  system. 
During  an  elastic  collision  the  energy,  T,  transferred  to  the  target  atom  is 

T  =  Tmsin2(G/2)  (3) 
where  Tm  is  the  maximum  energy  transfer 

Tm  =  Km2E^2'  (Id 
Finally,  the  connection  between  the  differential  cross  sections  in  the 
laboratory  system  and  in  the  center  of  mass  system  is 


6 \>j)   6(e)  P  +  8*+gyQ»e)* 


,  1+     cos  a  I 

The  differential  cross  section  for  scattering  in  a  given  direction 
B(O^)  is  defined  as  the  number  of  particles  scattered  into  a  solid  angle 
dS  per  unit  time  per  unit  incident  intensity  where  d£2  is  an  element  of 
solid  angle  in  the  direction  2.    For  azimuthal  symmetry  in  the  incident 
beam,  the  differential  cross  section  is  given  by 

d6-  2nbdb 
dS2  =  =  (6) 

where  b  is  the  impact  parameter.    The  element  of  the  solid  angle  d2  can  be 
written  as 

d2  =  2nsinQ  dO  (7 ) 

Substitution  of  equation  (7)  into  equation  (6)  gives 

6(Q)  =  _k__|db|  ,fl. 

5i.n9|d9|  (8) 

To  evaluate  equation  (8)  one  must  know  either  b(0)  or  9(b).    For  a  central 
potential  V  =  V(r),  the  scattering  angle  center  of  mass  system  and  the 
impact  parameter  are  related  by  the  formula 


6 


where  rQ  is  the  distance  of  closest  approach.     The  value  of  rQ  is  defined 
by  the  largest  positive  root  of  the  equation 

U        £.     /  '«     -  o    .  (10) 

Except  for  potentials  of  the  form  V(r)  -  c/rn,  equation  (9)  can  only  be 
solved  numerically.    In  addition,  difficulties  arise  if  the  potential  is 
attractive  (negative)  and  n  is  greater  than  two.    In  this  case  rQ  does  not 
exist  in  equation  (1)),  and,  in  general,  the  deflection  angle  is  then 
undefinable.    Since  all  potentials  discussed  in  this  thesis  are  repulsive, 
this  latter  problem  was  never  encountered.    The  integral  in  equation  (9) 
can  be  carried  out  in  closed  form  for  the  Coulomb  potential  V(r)  -  ZjZ?e2/r, 
giving  the  familiar  results: 

and 

[    J      V   4-E.  sin*  -f       /      '  (12) 

Following  Lindhard,15  another  potential  for  which  the  calculation  of  the 
scattering  angle  can  be  carried  out  in  closed  form  is 

where  s  is  a  -scaling  constant  and  a  is  a  constant  (dimension  of  length). 
The  differential  cross  section  is 

/n/   Tw'-'/n      -y  n-  •/»  1 
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T  is  the  energy  transfer  and  T    is  the  maximum  energy  transfer  defined  in 
equation  (I4). 

The  Coulomb  potential  includes  only  the  interaction  of  the  two  nuclei, 
and  the  screening  effects  of  the  electrons  cannot  adequately  be  described 
by  a  power  potential.    Other  potentials  have  been  proposed  to  examine  the 
screening  effect.    Section  III  will  discuss  some  of  these  potentials  and 
describe  the  technique  used  in  evaluating  equation  (10)  and  integrating 
equation  (9).    This  will  be  followed  by  the  presentation  of  the  numerical 
results  in  section  IV. 


III.    CALCULATIONAL  PROCEDURE 

There  have  been  many  potentials  proposed  which  estimate  the  screening 
effects  of  the  electrons.    To  investigate  the  effects  of  these  potentials 
on  the  classical  scattering  angle  and  the  differential  cross  section,  one 
must  numerically  solve  equation  (10)  for  rQ  and  numerically  evaluate  the 
integral  in  equation  (9).    We  discuss  below  the  different  forms  of  these 
potentials.    The  latter  part  of  this  section  will  deal  with  the  procedure 
used  in  calculating  the  classical  scattering  angle  and  the  differential 
cross  section. 

A .    Screened  Potentials 

To  approximate  the  screening  effects  of  the  electrons  in  interatomic 
interactions,  Bohr  proposed  the  potential 

V(r)  -  e-r/*  .  (13) 

The  screening  length,  a,  is 

where  aQ  is  the  radius  of  the  first  orbit  in  hydrogen.    The  screened  Bohr 
potential  decreases  too  rapidly1^  for  large  r  and  is  only  valid  in  the 
range  r/a0^l.    Attempts20  have  been  made  to  extend  the  range  of  validity 
by  revising  the.  value  of  the  screening  length,  a.    These  revised  values, 
however,  only  serve  to  extend  the  range  of  validity  by  a  very  small  amount. 

Thomas  and  Fermi9  developed  a  statistical  model  to  approximate  the 
screened  potential  for  an  atom.    The  form  of  the  Thomas-Fermi  potential  is 
V(r)  =  <j)(*).    The  screening  factor,  (j)(x),  is  determined  by  the 


? 


Thomas-Fermi  equation 

(15) 

<j>(o)=|     j  (j)(«>)  -  O        J       X  O 


The  riimensionless  variable,  x  =  ^ Ti  ovs —  >  allows  the  screening  function  to 
be  universal  for  all  atoms.    Since  equation  (15)  must  be  solved  numerically, 
|)(x)  is  expressed  in  the  form  of  a  numerical  table.10    For  the  interaction 
between  two  atoms,  the  Thomas-Fermi  potential  is  given  by 

0.8853ao  ii 
where  a-,v  =      ■  ■  -g-    —  .    Firsov-1"1-  proposed  a  similar  form  of  this 

potential, 

V(r;  -  <K^F)  (17) 

0.88^3a 

where  aF  =  ... .     r    Q      .    There  is  actually  very  little  difference  in 

equations  (16)  and  (17).  When  Z^/^  charigss,  the  ratio  aF/aTF  also  changes, 
but  this  variation  of  aF/aTF  is  less  than  +)4  per  cent,  and  usually  consider- 
ably less. -*-5 

Since  the  Thomas-Fermi  potential  is  not  expressed  in  a  closed  form, 
one  cannot  conveniently  use  this  potential  to  calculate  the  classical 
scattering  angle.    Moliere  proposed  an  analytical  curve  fit  to  the  Thomas - 
Fermi  potential.1-1    This  approximation 


(18) 


is  more  accurate  than  the  screened  Bohr  potential  set  at  large  r  and  has 
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often  been  used  in  channeling  calculations. 

Recently  much  work  has  bsen  devoted  to  approximate  analytical  solutions 
of  the  Thomas-Fermi  equation.    Csavinssky1^  has  proposed  the  trial  function 


where  a,  b,  c*  ,  and  &  are  parameters.    By  the  use  of  the  variational 
principle,  the  values  of  these  constants  (a  ■  0.711,  b  =  0.2889,        ■  0.175, 
and  /6  =  1.6625)  were  found  by  requiring  that  the  boundary  conditions 
(equation  (15))  be  satisfied,  a  +  b  =  1,  and  the  electron  density  be 
properly  normalized.    Csavinszky  states  that  this  potential  varies  from  the 
Thomas-Fermi  potential  less  than  one  per  cent  in  the  range  0  <  x  <  1.0  and 
less  than  eight  per  cent  in  the  range  1.0  <  x<2.0.     To  obtain  the  inter- 
action potential  between  two  atoms,  the  Firsov  approximation  (equation  (17 )) 
was  applied.    Firsov  limits  the  validity  of  this  approximation  to  the 
range  r/aQ  <  1.9. 

Another  trial  function, 


energy.    Their  criterion  for  a  better  potential  was  a  smaller  difference 
between  the  upper  and  lower  bounds.    They  also  checked  Csavinszky1 s  potential 
by  the  same  method  and  concluded  that  Roberts'  trial  function  was  a  better 
approximate  solution  to  the  Thomas-Fermi  equation. 


(19) 


(20) 


where  #   is  an  undetermined  constant,  was  suggested  by  Roberts. The 
optimum  value  of   $  =  1.7822  was  determined  by  Anderson,  Arthurs,  and 
RbbinsonAO  by  calculating  the  upper  and  lower  bounds  for  total  electron 


'See  for  example  C.  Erginsoy,  Phys.  Rev.  Letters  15,  ?60  (1965). 


Recently  Wedepohl1?  noted  that  the  Thorrias-Fermi  screening  function  can 
be  described  very  accurately  in  the  range  0.3  <  X  <  16  by  the  relation 

4)Cx>  =  Col7x)e-6-12x'4  (2D 

It  was  also  noted  that  this  form  of  <}>(x)  is  not  a  particular  solution  of  the 
Thomas-Fermi  equation  since  f)(0)  =  0  and  that  this  form  is  invalid  for  small 
r.    When  the  Firsov  approximation  is  applied,  the  interaction  potential  is 
given  by  the  following  relation 

V(r)  =  be'0"'17*  (22) 
where  B  =  1.55(10),J  Z1Z2ZQ  electron  volts,  c<  =  8.98ZQ  k~    t  and 
Z    =  4-  (Z^  +  zf*  )2.    The  difference  between  this  potential  and  the  Firsov 
potential  is  less  than  1.5  per  cent  in  the  range  0.3>  r/ap>l6. 

For  three  potentials  (Csavinszky,  Roberts,  and  Wedepohl),  the  values  of 
the  screening  functions  were  calculated  and  compared  in  table  I  to  the  exact 
Thomas-Fermi  screening  function  values  ."^    We  itfish  to  investigate  the 
behavior  of  the  classical  scattering  angle  and  the  differential  cross 
section  when  each  of  these  potentials  is  assumed.    The  numerical  methods 
used  in  these  calculations  will  be  discussed  in  the  second  part  of  this 
section. 

3.    Numerical  Procedure 

We  wish  to  examine  the  effects  of  the  five  screened  potentials 
(screened  Bohr,  Moliere,  Csavinszky,  Roberts,  and  .Jedepohl)  on  the  classical 
scattering  angls  and  the  differential  scattering  cross  section.    To  calculate 
the  classical  scattering  angle,  0,  and  the  differential  cross  section, , 
for  an  impact  parameter,  b,  three  numerical  evaluations  must  be  performed: 


TABLE  I 

Ratios  of  the  screening  functions  of  Csavinszky,  <£>c, 
Roberts,   <p^  ,  and  Wedepohl,   <£w  ,  to  the  Thomas-Fermi 
Screening  function,   tyx)  ♦ 
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2.60 
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1-.18U10 

0:99667 

2  ;80 

0-.169880 

1.13129 

I.I880I1 

0.99776 

3.00 

0.156630 

1-.1U035 

1.19099 

0.99897 

T.20 

0.iLiL820 

1.11:720 

r.19293 

1; 00019 

3".  h0 

0.13h250 

1-.15180 

1.19392 

1.00139 

t.6o 

0.12U7UO 

r.i5h28 

1.19U1U 

1-.  00265 

3 '.80 

0.116170 

1.15I150 

1.193U8 

1:00380 

U-.oo 

0.108)400 

1.15276 

1.19222 

1-.0050U 

L.60 

O.OP90R5 

1.13558 

1-.18U17 

1.00816 

5 -.00 

0.078808 

l.H55!i 

1.17595 

1.00992 

6.00 

0.059H23 

1.0U216 

1.1U7U9 

1.01296 

7-00 

O.OL16098 

0.9)4660 

1.11056 

1.01396 

8.00 

0.036$8? 

0.81.0U5 

1.06799 

1.01302 

9.00 

0.029591 

0.7^228 

1.0218U 

1.01029 

10.00 

0.02h3iLi 

0.62802 

0.97  372 

1.00600 

TABLE  I  (continued) 


4>Cx)        4>c/^       ^  4V(j) 


0.010805  0.2U558  0-.73527  O.967I6 
o.oo?h7li       o-.o?"»07       0-.38U87  0.85003 

0.0006? 2         0.00002         0.072^9  0.$6755 


Ill 


1)  solving  equation  (10)  for  the  distance  of  closest  approach,  r  , 

2)  evaluating  the  integral  in  equation  (9),  and  3)  evaluating  db/dO.  The 
methods  used  to  make  these  three  evsluations  will  now  be  discussed. 

1.    Evaluation  of  rQ.    The  distance  of  closest  approach,  r0,  is 
defined  by  the  equation 


For  a  repulsive  potential  the  left-hand  side  of  this  equation  is  a 
monotonically  increasing  function  of  rQ  when  l>V(rQ)/E.    Except  for 
collisions  which  cause  large  deflections,  the  values  of  rQ  and  b  vary  only 
by  a  smell  amount.    The  value  of  ro  is  found  by  locating  the  point  of 
intersection  of  the  curve  f(x)  =  (1  -  ( V(x)/E)  J^x  and  the  straight  line 
g(x)  =  b.    The  operation  of  a  mesh  used  to  locate  this  point  of  intersection 
is  described  by  the  following  steps.    1)  First  starting  with  x  =  b  and 
increasing  x,  if  necessary,  find  a  value  of  x  (x  =  s)  where  (1  -  (V(s)/E)V2 
is  real.    This  s  is  the  starting  value  for  the  mesh.     2)  The  value  of  f(s) 
is  checked  and  found  less  than  (greater  than)  b.    Then  s  is  repeatedly 
increased  (decreased)  by  the  mesh  size  until  f(s)  is  greater  than  (less  than) 
b.    At  this  point  we  know  the  value  of  rQ  accurate  to  the  value  of  the  mesh 
size,    h)  The  mesh  size  is  reduced  and  steps  2),  3),  and  i|)  are  repeated 
until  rc  is  known  to  the  accuracy  desired.    If  a  mesh  point  ever  coincides 
with  the  point  of  intersection  of  the  two  curves,  the  value  of  rQ  is  then 
known  and  no  further  calculations  are  needed.     In  the  present  calculations, 
the  mesh  size  was  varied  from  0.01aQ  to  (10~7)aQ 

2.    Evaluation  of  0.    The  classical  scattering  angle,  9,  for  a 
collision  between  two  atoms  is  given  by 

e  -  U  -  2 b  /  °°r  11 1  -  r *  -  bS> &        '  (9a) 


15 


One  difficulty  arises  in  a  numerical  evaluation  of  this  equation.  Since 
there  exists  a  singularity  at  r  =  rQ  in  the  integrand,  the  numerical 
integration  can  run  into  trouble  at  this  point.    Recently  Smith"  suggested 
the  use  of  the  Gauss-Mehler  quadrature  to  avoid  this  difficulty.    If  we 
substitute  x  =  rQ/r  into  equation  (9a),  it  becomes 

e(b)-T-2k  (*  £<&A* —  (23) 


where 


frx)  = 


Since  the  singularity  behaves  as  (l  -  x^)"""^  at  x  =  1,  it  is  possible  to 
remove  this  singularity  by  the  use  of  a  special  Gaussian  quadrature 
f ormula22  derived  by  Mehler  for  integrals  which  have  a  weighting  function 
of  (1  -  u2)-l/2.    This  quadrature  formula  is 

f-  f(u)(l-U^Ju=f  If  (cos  1T)  • 

Since  the  potential  is  spherically  symmetric,  f(x)  =  f(-x),  and  equation  (23) 
can  be  written  as 

This  can  be  evaluated  by  using  the  Gauss-Mehler  formula, 


J-  I 

If  n  is  even. 


and  the  angle  0  can  be  approximated  by  the  series 
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This  method  is  very  effective  and  is  well  suited  for  a  computer.  Smith 
calculated  the  deflection  angles  for  collisions  of  two  hydrogen  atoms, 
assuming  an  interatomic  potential  of  V(r)  =  d/r^,  and  obtained  an  accuracy 
of  three  significant  figures  with  only  two  integration  points  (n  =  U). 
Similar  accuracy  is  found  for  other  repulsive  potentials  provided  that  f(r) 
does  not  come  close  to  zero  for  any  value  of  r  >rQ.    This  condition  is 
satisfied  for  all  five  potentials  discussed  in  the  first  part  of  this 
section.    Following  the  suggestion  of  Tang  and  Karplus,^3  n  was  set  equal 
to  thirty  in  the  present  calculations. 

^.    Evaluation  of  & (Q) .    The  differential  cross  section,  6"(8),  is 
determined  by  the  equation 

"iiiri"  pfel  (8a) 

Siace  Q(b)  cannot  be  expressed  in  a  closed  form,  dQ/db  must  be  numerically- 
evaluated.    This  was  accomplished  by  a  four  point  central  difference 
f ormula 


(21;a) 


where  h  =  for  all  i  (i.e.,  all  points  must  be  equidistantly  spaced) 

and  the  error  is  approximately  IrV    The  derivatives  at  the  first  two  points 
were  found  by  a  forward  difference  formula 

and  the  derivatives  at  the  last  two  points  were  found  by  a  backward 
difference  formula 

^L=IHy;-3  +  fy^-3y^^y;]+  &w>  (2Uc) 
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rf-(Q)  for  various  values  of  b,  ranging  from  C.001^ao  to  2.5aQ,  were 
calculated  for  each  of  the  five  potentials  (screened  Bohr,  Moliere, 
Csavinszky,  Roberts,  and  VJedepohl).    The  impact  parameters  were  separated 
into  three  groups:     one  corresponding  to  small  b  (0.0015ao  <  b  <  0.0)48ao) 
with  h  =  O.OOl^a    (32  points),  another  corresponding  to  central  b 
(0.05a0<  b  <  0.6ao)  with  h  =  O.25a0  (19  points),  and  the  third  corresponding 
to  large  b  (0.?aQ^  b  <  2.$a0)  with  h  =  0.1ao  (21  points).    At  each  value 
of  b,  a  value  of  db/dQ  can  be  found  by  equations  (2U)  and  the  differential 
cross  section  for  this  value  of  b  can  then  be  calculated  by  formula  (8). 

In  this  section  we  examined  the  potentials  proposed  by  Bohr,  Moliere, 
Csavinszky,  Roberts,  and  VJedepohl.    In  the  latter  part  of  this  section  the 
numerical  methods  used  to  calculate  the  distance  of  closest  approach,  the 
classical  scattering  angle,  and  the  differential  cross  section  were 
discussed.    In  section  IV,  we  will  show  the  reliability  of  these  numerical 
methods  by  comparing  present  calculations  to  the  calculations  made  by  Everhart 
assuming  the  screened  Bohr  potential.    Results  of  calculations  for  each  of 
the  five  potentials  will  then  be  presented  followed  by  a  discussion. 


IV .    NUMERICAL  RESULTS 

The  classical  scattering  angles  and  the  classical  scattering  cross 
sections  were  numerically  evaluated  by  the  methods  discussed  in  section  III 
for  five  screened  potentials  (Bohr,  Moliere,  Csavinszky,  Roberts,  and 
Wedepohl).    To  show  the  validity  of  the  numerical  methods,  present  cal- 
culations for  the  screened  Bohr  potential  were  compared  to  similar  cal- 
culations of  Everhart,  Stone,  and  Carbone.3    Tne  numerical  results  of  the 
scattering  and  cross  section  calculations  for  each  potential  will  then  be 
presented. 

A.    Comparison  with  Previous  Calculations  (Bohr  potential).  Everhart, 
Stone,  and  Carbone8  presented  in  a  tabular  form  their  evaluations  of  the 
distance  of  closest  approach,  rQ,  the  impact  parameter,  b,  and  the 
differential  cross  section, 6  (9) ,  for  different  values  of  scattering  angles, 
Q.    These  calculations,  made  for  various  values  of  f/a  where  a  is  the 
screening  length  and    f  is  the  collision  diameter,  were  made  universal  by 
expressing  the  results  in  dimensionless  quantities:  rQ/p  ,  b/p  ,  and 
6(Q)/p  2.    The  collision  diameter,  P  =  ZXZ2§  /E  where  E  is  the  energy  in 
the  center  of  mass  system,  is  the  distance  of  closest  approach  for  a  head-on 
collision  (no  screening).    To  compare  the  calculations,  a  value  was  assigned 
to       and  Z2  (thus  fixing  a).    For  a  value  of  f /a,  the  values  of  b/p 
were  used  to  specify  E  and  b.    These  values  of  energy  and  impact  parameter 
were  used  to  calculate  r    and  9,  and  then  the  values  of  Q  and  rQ/  f>  were 
compared  to  those  of  Everhart,  Stone,  and  Carbone  for  all  values  of  P  /a 
(except  (r /a  =  0)  which  corresponds  to  Rutherford  scattering).  Their 
results  varied  from  the  present  results  by  an  average  of  0.5  per  cent. 
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Two  of  these  comparisons  are  presented  in  table  II.    This  comparison  was 
carried  out  for  four  different  sets  of  values  for  Z±  and  ?2  and  it  was  found 
that  the  calculations  varied  by  a  maximum  of  0.06  per  cent.    The  differential 
cross  sections  tf"(9)/p    ,  were  also  compared,  and  again  closs  agreement 
(variational  per  cent)  was  found.    In  addition,  calculations  were  made  to 
extend  the  tables  to  include  small  angle  scattering.    These  are  presented 
in  tables  III,  IV,  V,  and  VI.    Now  that  the  reliability  of  this  method  is 
established,  we  are  ready  to  investigate  the  effects  of  the  different 
potentials  on  the  scattering  angle  and  the  cross  section. 

B*    Results  of  Calculations  Comparing  Potentials.    The  scattering  angle,  9, 
was  calculated  for  each  cf  the  five  potentials  (Bohr,  Koliere,  Csavinszky, 
Roberts,  and.  Wedepohl)  for  several  values  of  energy,  Z-^  and  Zg.    Results  of 
two  cases  (E  =  21. I4.  kev  and  E  =  [}2.9  kev)  are  presented  in  tables  VII  and 
VIII  for  Z-j_  =  8  and  Z2  =  6.     To  aid  comparison,  for  each  value  of  b,  9R 
calculated  using  the  Roberts  potential  is  shown  along  with  the  ratios  of 
the  angles  to  9R  calculated  using  the  other  potentials. 

The  differential  cross  section  was  also  calculated  in  each  of  the 
above  cases.    For  small  angles,  the  calculated  differential  cross  section 
corresponding  to  the  above  cases  are  given  in  tables  IX  and  X.    To  represent 
£(Q)  in  an  universal  form,  Lindhard,  Nielsen,  and  Scharff^  introduced 
the  parameter  t  given  by 

i  =  e-  sin1  .a    ;      £  =  a'/p  (25) 

where  a' =  0.3853a  (a  is  the  screening  length  defined  in  equation  (1U)). 
The  cross  section  then  was  written  in  the  form 


TABLE  II 


Comparison  of  the  present  calculations  of  0  and  r0/ f     =  A  with  those 
of  Everhart  et  al.°    The  ratios  are  given  in  columns  h  arid  6.    b  and 
r0  are  in  units  of  aQ  =  0.53(10"°)  cm  and  9  is  in  degrees. 


P /a  =  $-0,  Energy  =  lhh'b  ev  (center  of  mass  system) 


b 

ro 

9 

9E/9 

A=r0/r- 

ratio 

2.O667676 

2.0782 

1.81060 

0.99U15 

0.91906 

0.99993 

1.78UH3U 

1.8051 

3.57663 

1.00653 

0.79829 

0. 9996U 

1. 6I67 821 

l-.6b.67 

5  •  3h0!40 

1.01116 

0.72823 

0.99831 

I.I.969367 

1.5^3 

7.10329 

1.01*61 

0.67398 

1.00002 

1.^27 3UU2 

l.*8l7 

10.596*6 

1.01922 

0.61103 

0.99831 

1.2007152 

1.2707 

1U.2321U 

r.  01179 

0.56196 

0.99829 

1.1012208 

1.1863 

17.8Q913 

1.0056)4 

0.52h6h 

1.00069 

0-929*671 

I.OU76 

26.hU 07 6 

1.02115 

0.U6327 

1.0037)4 

0'.79595b3 

0.9U72 

*5 -.61025 

1.0109U 

O.U1891 

1.00261 

0.60*7b95 

0.8172 

5U-.27U2U 

0.99)495 

O.36IUO 

0.99612 

0.U725978 

0.7U09 

71.98321 

1.00023 

0.°276h 

I.OOUII4 

0.366*199 

0.688b 

90.ID485 

0.99873 

O.3OUIJ4 

1.00185 

0.1650701 

0.6191 

135.00760 

0.9999b 

0.27381 

0.9970U 

0.0002261 

0.6000 

179.93636 

1.00035 

0.2653)4 

0.99870 

r/a  =  0.2, 

Energy  =  3610 

ev  ( center  of 

mass  system) 

b 

ro 

9 

9£/9 

A=r0/r 

ratio 

0.3601733 

0.8669 

1.78958 

1.00582 

9-53383 

I.OOO6I4 

0.63221,09 

0.6U32 

3.5732I4 

1.007)1 9 

7.11162 

0.99977 

O.5083253 

0.5228 

5.39*30 

1.0012)4 

5-77959 

0.9983b 

O.U2963UI4 

'  0-.LU68 

7.162*8 

1.00525 

b. 93981 

1.0000b 

0.*283311 

0.3)499 

10.7).i  80li 

1.00U83 

3.86800 

1.00052 

0.26U1121 

0.2891 

1U.U3917 

0.99729 

3.19577 

0.99819 

0.220696h 

0.2U8U 

18.10185 

0.99b 37 

2.7b57b 

0.99791 

0.1$6h773 

0.1891 

26.97016 

1.00111 

2.09020 

0.99990 

0.1557 

*5.95565 

1.00115 

1.72113 

1.00515 

0.0732*87 

0.1202 

5b. 0*629 

0.99933 

1.32937 

1.00799 

O.O556263 

O.IO2I4 

71.97532 

1.0003)4 

1.1*167 

0.99852 

TABLE  II  (continued) 
p/a  =  0.2,  Energy  =  3610  ev  (center  mass  system) 


9g/9  A=r0/p  ratio 


0.0UO7O22        0.0919  89'.  91896        1.00068        1.01558  1.00U35 

0-.01691U0  0.0795  1^5.09512  0.99930  O.8786U  0.99927 
0.0000090        0.076U        17 9. 97 170        1.0001U        0.8Uh58  1.00050 
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TABLE  III 


6(®)/r  2  f°r  different  values  of  impact  parameter,  b,  for 
the  screened  Bohr  potential.     p  =  Z^Zo^/e  •    9  is  in 
degrees.    Hie  nxanber  in  parentheses  after  an  entry  is  the 
power  cf  10  to  be  multiplied  by  it. 


p/a  -  0.1 


67.5552  67.5558 

6b/.  85  30  6U.8537 

62.1508  62.1518 

59.Ub.85  59.UU99 

56.7U63  56-7U81 

5U.0UU1  51-OU6U 

51.3U19  51-3UU9 

LB-.6397  U8.6U36 

U5-9375  1*5  91*26 

U3-2353  U3-2U19 

h0.5~3l  U0.5U18 

37 -.8309  37-8U23 

35  ■•1287  35 '.11*36 

32-.U265  32-.J.U60 

29.7  2U3  29.71-98 

27 -.0221  27.0555 

2U-.3199  2U.3637 

21.6177  21.6750 

18 -.9155  18.9905 

16-.2133  16.311U 

13.5m  13.6395 


0 

6  {Q)/p  c 

0.003U 

0.11*21(12; 

0.00U6 

0.1251(12) 

0.0060 

0.5239U-U 

0.0082 

0.268U(11) 

o.muo(ii) 

0.0109 

0.O1U7 

0.7385(10) 

0.0198 

0.3853(10) 

0.0268 

0.1975(10) 

0.036U 

0.1012(10) 

0-0U93 

0.5138(  9) 

O.0672 

0.2588(  9) 

0.0915 

0.1280(  9) 

O.1255 

0.6296(  8) 

0-.1719 

0.3069(  8) 

O.2369 

0.1UU9(  8) 

0.3287 

0.6756(  7) 

0.h579 

0.3063(  7) 

0.6UU3 

0.1336(  7) 

0.9161 

0.5596 (  6) 

1.3210 

0.2235(  6) 

1.9b ho 

0.8109(  5) 

f/a  =  0.2 


b/p  r0/p  0  6(Q)/p2 


36.h798  36.1*801  0.0037  0.3703(11) 

35;1287  35.1291  0.0050  0.2261(11) 

33.7776  33.7782  O.OO67  0.1110(11) 

32.U265  32.1*272  0.0091  0.5929(10) 

?1.075h  31.0761*  0.0121  0.3110(10) 

29.721*3  29.7256  0.0165  0.1662(10) 

28.3732  28.37U9  0.0219  0.8883(  9) 


TABLE  III  (continued) 


/a  =  0.2 


t>/r 


27  '.02  21 
25.6709 
214.3199 
22.9688 

21.6177 
20.2666 

18.9155 

17  -56U3 
16.2132 
lh.862l 
13.5110 
12.1599 
10.8088 

9-U577 


rj., 


r 


27.02b? 
25.67^9 
2U-3237 
22; 97 38 
21.621*3 
20.2752 
18.9268 
17.5792 
16.2327 
IU.8876 

13'^UhU 
12.2036 
10.8659 
9.5323 


Q 

2 

0.0296 

o.)i5hi( 

9) 

o.okoo 

0.2L01( 

9) 

0.05^8 

o.i2U5( 

9) 

0.0728 

0.6330( 

8) 

0.0987 

0.3211( 

8) 

0.13U3 

0.1616( 

8) 

0.1830 

0.8055( 

7) 

0.2502 

0.3968( 

7) 

0.3U30 

0.1926( 

7) 

O.U72U 

0.9159( 

6) 

0.65U1 

O.U270( 

6) 

0.9113 

0.19U1( 

6) 

1.2793 

0.8608( 

5) 

1.8D.I 

0.3632( 

5) 

TABLE  TV 


6(Q)/r  ?  for  different  values  of  impact,  parameter ,  b,  for 
the  screened  Bohr  potential,     p  =  2]_22e  /£  •    0  is  in 
deerees.    The  number  in  parentheses  after  an  entry  is  the 
power  of  10  to  be  multiplied  by  it. 


T/a  =0.5 

b/jp  r0/p 

15  -.67  2  8  15-6730 

15 -.1321*  15.1326 

DV.5919  lk.5923 

m-.0515  l!i.0519 

13.5110  13.5116 

12.9706  12.97U4 

12.U^02  12.U312 

11.8897  11.8910 

11.3U93  11.3510 

10.8088  10.8111 

10. 268)4  10.2713 
9. 7280  9-7318 
9.1875  9.1926 
8-.6U71  8.6537 
8.1066  8.1153 
7-5662  7.5775 
7.0257  7.0UG6 
6.U853  6.50U7 
5-.5UU9  5.9702 
S-hohh  5-U375 


,36UO  h.9072 


9 

6(Q)/r 

2 

0.0053 

0.3930(10) 

0.0070 

0.1790(10) 

0.0095 

0.9592( 

9) 

0.0127 

0.52)47  ( 

9) 

0.0170 

0.2BU0( 

9) 

0.0228 

0.1537 ( 

9) 

0.03O14 

0.8212 ( 

8) 

0.0U08 

O.U286( 

8) 

0.05U9 

0.2259( 

8) 

0.0738 

0.1182( 

8) 

0.0996 

0.6131( 

7) 

0.13UU 

0.3180( 

7) 

0.1817 

0.1629( 

7) 

0.2U6U 

0.8290( 

6) 

0.33U5 

O.U190( 

6) 

O.U553 

0.209U( 

6) 

0.6213 

0.103U( 

6) 

0.8508 

0.5035 ( 

5) 

1.1689 

0.2105  ( 

5) 

1.613U 

0.11U3( 

5) 

2.2391 

0.52U7( 

k) 

/a  •=  1.0 

b/p  r0/r  9  6{l)/p 


8-. 3766  8.3770  0.0059  0.7281(  9) 

8.1066  8.^068  0.0079  0.h099(  9) 

7'.P^6U  7.8366  0.0105  0.21!41(  9) 

7.5662  7-566)4  0.01U1  0.115K  9) 

7.2960  7.2963  0.0189  0.6271(  8) 

7.0257  7-0262  0.0252  0.^3k6(  8) 


TABLE  IV  (continued) 


f/a  - 


1.0 


6.7555 
6.1*653 
6.2151 
5 -.91*1*9 

5 -.1*01*1* 

5.131*2 

1*-.  861*0 
h.5938 
U-3235 
h.0533 
V.7831 

3.5129 
3.21*27 
2.972U 


-0/ 


r 


6. 7561 

6.1*861 

6.2163 

5-9U62 

5.6763 

5-U067  . 

5.1371 

h.8678 

/V-5988 

b.3^01 

!*.06l9 

■*:79hh 

3.5276 

^.2619 

2.9975 


0.0339 
0.01*51* 
0.0608 
0.0816 
0.1098 
0;  11*76 
0.1987 
0.2681 
0.?626 
0:1*907 
0.6650 
0.9036 
1.2305 
1.6797 
2.2997 


0;l77l*(  8) 

0.95U5 (  7) 

0.5105C  7) 

0.2680(  7) 

0: 11*10  (  7) 
0'.7l*22( 
0.3867( 

0:1995( 
0.1026( 

o-.5258(  . 
0.2661*(  5) 

0.1336f  5) 

0. 661*5  ( 
0.^28)*( 

0.1589 ( 


6) 
6) 
6) 
6) 
5) 


I*) 

U) 
I*) 


TABLE  V 


6{Q)/pd  for  different  values  of  impact  parameter,  b, 
for  the  screened  Bohr  potential,    p  =  2tZ2 &  /e  •    9  is 
in  degrees.    The  number  in  parentheses  after  an  entry  is 
the  power  cf  10  to  be  multiplied  by  it. 


p/a  =  2.0 


b/p 

V/"1 

0 

/  t r,\  I „  2 

6  W/p 

]y.Li586 

Li.U587 

0.0068 

U.3235 

I4.3236 

0.0089 

0.  foyo{ 0) 

h.loou 

Li  .loop 

J  •  1X1X7  7  V  O  / 

b-0533 

h.0535 

0.0159 

0.2h82(8) 

y.  9132 

r.918h 

0.0212 

0.13U7(8) 

3.7  831 

3.7833 

0.0232 

0.7256(7) 
0.3870(7) 

3-.6U30 

^.6U83 

0.0378 

3.5129 

3.5133' 

0.0506 

0.2096(7) 

3.3778 

3;  37  83 

O.O676 

0.1123(7) 

3 -.2)426 

3.2U3U 

0.0907 

0.5972(6) 

3 -.1075 

3.1085 

0.1215 

0.3192(6) 

2-.972U 

2.9737 

0.1629 

0.1690(6) 

2.8373 

2 ;8390 

0.2188 

0.8880(5) 

2.7022 

2.70l4h 

0.29UU 

o.U665(5) 

2.5671 

2.5700 

0.3961 

0.2U!jO(5) 

2.l|320 

2.U358 

0.5337 

0.1263(5) 

2.2969 

2.3019 

0.7199 

o.6558(U) 

2.1618 

2.1683 

0.9723 

0.3367(U) 

2.0267 

2.0352 

1.3151 

0.1719(U) 

1.8915 

1.9027 

I.78OU 

0.8769(3) 

1.7561* 

1.7710 

2.hl38 

0.Uh08(3) 

P/a  =  5-0 


1-.9U56  1.9h56  0.0071  0.2336(8) 

1.3915  1.3916  0.0095  0.12hU(8) 

1.3375  1.8376  0.0126  0.7290(7) 

r.7835  1.7835  O.OI67  0.39h7(7) 

1-.729U  1.7295  0.022b  0.2109(7) 

1.675)1  1-67  5  5  0.029  3  0.1171(7) 

1.6213  1.6215  0.0^96  0.6!i09(6) 


TABLE  V  (continued) 


p/a  =  5-0 


ro/r 

9 

1.5673 

1.5675 

0.0529 

0.3U37(6) 

1.5132 

1.5135 

0.0707 

0.185U(6) 

1.U592 

1-U595 

0.09UU 

0.1001(6) 

1.U051 

1.UQ56 

0.1262 

0.5335(5) 

1-.3511 

1.3517 

0.1688 

0.2898(5) 

r.2971 

1.2978 

O.2257 

0.15U8(5) 

1-.2U30 

1.2UU0 

0.3025 

0.82U3(U) 

1.1390 

1.1903 

O.I4O52 

O.U386(U) 

1-.13U9 

1.1366 

O.Sh3k 

0.2323(h) 

1.0809 

1.0831 

0.7290 

0.1227(h) 

1.0268 

1.0297 

0.9786 

0.6h62(3) 

0.9728 

0.9766 

1.3138 

0.3393(3) 

0-.9188 

0.9237 

1.76UU 

0.1733(3) 

O.86U7 

0.8711 ' 

2.3689 

0.9332(2) 

TABLE  VI 


^(Q)/p     for  different  values  of  impact  parameter,  b, 
for  the  screened  Bohr  potential.    P*  lt2xe.A/lZ.      •  9 
is  in  decrees.    The  number  in  parentheses  after  an 
entry  is  the-  power  of  10  to  be  multiplied  by  it. 


p/&  *  10.0 


b/r 

9 

1.0268 

1.0269 

0.0081 

O.U989(7) 

0.9998 

0.9998 

0.0107 

0.2700(7) 

0.9728 

0.9728 

0.01U2 

0.11*90(7) 

0.9158 

0.9)458 

0.0190 

0.8152(6) 

0.9188 

0.9188 

0.0252 

0.h529(6) 

0.8917 

0.8918 

0.0335 

0;  2)466(6) 

0.86U7 

0.861*8 

0.0UU7 

0.13U0(6) 

0.S377 

0.8378  ' 

0.059U 

0.73U0(5) 

0.8107 

0.8103 

0.0793 

0.3963(5) 

O.7836 

0.7838 

0.1058 

0.2155(5) 

0.7566 

0.7569 

0.1U11 

0.1172(5) 

0.7296 

O.7299 

0:1883 

0.6315(10 

0;7Q26 

0.7030 

0.2515 

0.3U01(U) 

0.6756 

O.676I 

0.3361 

0.1828(h) 

0.6U85 

0.6U93 

0-.!iU9U 

0.9807(3) 

0.6215 

0-.6225 

0,6009 

0.5258(3) 

0-5958 

0.8035 

0.2811(3) 

06&6 

0.5692 

1.07U6 

0.1501(3) 

0.8009(2) 

0-.5U0U 

0.5U26 

1-.U36U 

0.513a 

0.5163 

1.9196 

O.U283(2) 

O.U86U 

0.L901 

2.5628 

0.2288(2) 

For  Zi  -  8,  Z2  -  6,  and  E  =  1*2, 800  ev  (ELab  -  100  kev)  cal- 
culated scattering  angle  with  Roberts  potential,  Qo.  The 
impact  parameter,  b,  is  in  units  of  aQ  =  0.53(lO-°)  cm. 
Columns  3-6  represent  the  ratios  of  the  scattering  angles 
calculated  with  the  Bohr  potential  ( B) ,  the  Csavinszky 
potential  (C),  the  Moliere  potential  (K),  and  the  Wedepohl 
potential  (W)  to  the  corresponding  value  of  0j>.    0  is  in 
degrees.    The  number  in  parentheses  after  an  entry  is  the 
power  of  10  to  be  multiplied  by  it. 


b 

9R 

&  R 

C'  R 

On/On 

&  j/0o 
w'  rt 

2.5000 

0-U97K 

-1) 

0.0568 

0.8576 



1.0 8U9 

O.9OI46 

2 .  U000 

0.5616( 

-1) 

O.O67I 

0 . 87  80 

1.0821 

0 . 39U0 

2 ; 3000 

0.63 33( 

-1) 

0.07  9)i 

0.900U 

1.0817 

0.8891 

2. 2000 

0.7 200 ( 

-1) 

0.0922 

0.9160 

1 . 07 11 

O.878O 

2 . 1000 

0.8183( 

-1) 

0.1103 

O.933U 

l;066h 

0.87  02 

2 .0000 

0-9357(-D 

0.1296 

0.9)i90 

i.0561 

0.8601 

1 • 9000 

0.1071 

0 . 96  37 

l  .0^0)1 

0.856U 

1 . 8000 

0.1235 

0.1800 

0 . 97  K  3 

I.O378 

O.8502 

1.7000 

0.1U30 

0.2091 

0.98141 

1.0256 

0.81,37 

0.1665 

0.2)i)i7 

0.989^ 

1.0119 

0.837)1 

l'.cJOOO 

0-.19U7 

0.283)4 

0.99)i3 

0.9982 

0-83)4.6 

l.liOOO 

0.2292 

0.3295 

0.99U8 

0.9860 

0.8322 

1.3000 

0.2718 

0-.3798 

0.9936 

0-97U2 

0.8301 

1.2000 

0.3252 

O.U352 

O.9876 

O.962I4 

0.8300 

1.1000 

0.3923 

0-U973 

0.9792 

0.952h 

0.831U 

1.0000 

0-U791 

0.56U6 

0.9672 

0.9Uh6 

0.8339 

0.9000 

0.5925 

O.6369 

0.9526 

0.9U05 

0.8385 

0.8000 

0.7hU5 

0.7130 

0.9365 

0.9396 

0.8U47 

0.7000 

0.9539 

0.7917 

0.9198 

0.9U26 

0.85U1 

0-.6000 

0.1252( 

1) 

O.8703 

O.9063 

0.9510 

0.8676 

0.5000 

0.1699( 

l) 

0.9UU8 

0.8955 

0.9613 

0.8337 

0.L750 

0.18UU( 

1) 

0.9623 

0.8951 

0.9650 

O.8887 

0.L500 

0.2010( 

1) 

0.979U 

0.8939 

0.9678 

0.893U 

0.it250 

0.2197 ( 

1) 

0.9957 

0.89U7 

0.971U 

0.8989 

0.1-000 

0. 21*11  ( 

1) 

1.0108 

0.8957 

0.97U9 

0.9050 

0-.37S0 

0.2657( 

1) 

1.02U9 

0.8969 

0.97 81 

0.9107 

0-.3500 

0.29h2( 

1) 

1.0381 

0.9000 

0.9315 

0.9169 

0.3250 

0.3276( 

V. 

1.0501 

0.9037 

0.98U3 

0.9236 

0.3000 

0.3671( 

1) 

1.0603 

0.9087 

0 . 9869 

0.9302 

0.2750 

0.iaU2( 

l) 

1.0697 

0.91U7 

0.9B97 

0.9378 

0.2500 

0-h7l7( 

1 ) 

I.0767 

0.9219 

0.9916 

0.9)4li6 

0.2250 

0-5U26( 

1) 

1.0819 

0.930^ 

0.9933 

0.9513 

0.2000 

0.6320( 

I) 

1.0853 

0.9397 

0.99h8 

0.9580 

0.1750 

0.7U79( 

1) 

1.0360 

0.9503 

0.9961 

0.9639 

TABLE  VII  (continued) 




b 

Od 

— —  

Qr/0t? 



©h/Qd 

S-.t/Qr 
aiU  "it 

0.1500 

 ,  

0.9033( 

I) 



1.03U7 

0.9618 

0.9975 

0.9687 

0.1250 

0.1122( 

2) 

1.0800 

0.97  29 

0.998^ 

0."70k 

O'.IOOO 

0.1h50( 

2) 

1.0728 

0.93h^ 

1.0002 

0.Q679 

0.0750 

0.199l( 

2) 

1.062U 

0.99li5 

1.0018 

0.9^7  3 

0;Of>00 

2) 

1.0)18^ 

1.0016 

1'.0035 

0.9317 

0.  oil  80 

0.317  3( 

2) 

1.0022 

1:0036 

0:9288 

0'.  0)1 6s 

0.327^ ( 

2 ) 

1.0)  16)1 

1.0025 

l.OOlj.0 

0.9266 

0.0[j50 

0.3385( 

2  ; 

1.0)i53 

1.0026 

1. 0038 

0.92)il 

0'.3500( 

2) 

1.0,i)i3 

1.0029 

I.O039 

0 . 9217 

0.0l|20 

0.36  2)4  ( 

2) 

1.0li32 

1.0033 

l .  00)  1 1 

0.9T  90 

0 . 0h05 

0. 37  tt( 

2) 

1.0li23 

1.003U 

1  '.OOlil 

0.9l6)i 

0.0390 

0: 389^( 

2) 

1  •.  0)1  iii 

1.0038 

1  00)1? 

-1-  •  WU/4  c 

0  91  ?7 

O'JiOlili  f 

2) 

1  '.0)i03 

1  ■  nn)i)« 

n  9ii n 

0.0360 

O.U206( 

2 ) 

1.039? 

l'.00[iO 

1  ■  nolil 

X  •  \_'  \J  L\.  JL 

0:  9079 

n.03'is' 

0.)iT7R( 

2) 

1.0380 

1 . 00)  1 2 

n'  on^i 

0.0^30 

0.)i^6^( 

2 ) 

1.0368 

1.00)42 

1  . 00)  1  ? 

0.9018 

0.0315 

0:h765( 

2) 

1'.0358 

1.00U5 

1.00h3 

0.8° 89 

0.0300 

0.h982( 

2) 

1.031,6 

1;00UU 

1  •  OOli 

w  •  U  7j7  O 

0.0285 

0.^217 ( 

2 ) 

1.0^5 

l.OOli'i 

n  R9?9 

C.O27O 

2) 

1.0322 

l .  00)  1 6 

1  '.OOlil. 

0:8897 

0.02^ 

0.^7^2( 

2 ) 

1.0309 

1.00)t6 

1 .OOh^ 

0.^867 

0.02h0 

0 . 6057 ( 

2) 

1.0297 

1 .  OOlj.6 

1  on)  1 1. 

0  RRlil 

0.0225 

0.6392( 

2) 

1.0282 

I.OOI4I 

0.881^; 

0.0210 

0.676o( 

2) 

1.0268 

1  .00)i^ 

n  8791 

0 .019^ 

0.7 16^ 

2) 

1.  OOlil 

1  on^9 

n  877  ^ 

0.0180 

0.7612( 

2) 

1.0238 

l.OOhl 

1.0038 

O.8760 

0.0165 

0.8108( 

2) 

1.0223 

1.0038 

1.0036 

0.8755 

0.0150 

0.8656( 

2) 

1.0207 

1.0036 

1.0035 

O.8762 

0.0135 

0.9266( 

2) 

1.0188 

1.0033 

1.0032 

O.8780 

0.0120 

0.99U3( 

2) 

1.0170 

1.0031 

1.0030 

0.8817 

0.0105 

O.IO69 ( 

3) 

1.0150 

1.0028 

1.0027 

0.8875 

0:0090 

0.1152( 

3) 

1.0130 

I.OO2I4 

1-.002U 

0.8957 

0.0075 

0.12U3( 

3) 

1.0107 

1.0020 

1.0020 

O.9069 

0.0060 

0.13)i2( 

3) 

1.0085 

r.0016 

1.0016 

O.9213 

0.00)45 

0.1LU8( 

3) 

1.0063 

1.0012 

1.0012 

0.9388 

o.oo^o 

0.1$62( 

3) 

1.00)40 

1.0008 

1.0008 

0.9587 

0.0015 

0.1680( 

3) 

1.0019 

1.000U 

1.000b 

0.9797 

TABLE  VIII 


For  Z1  =  8,  Z2  -  6,  and  E  =  21,1*1*0  ev  (Exab  =  50  kev)  cal- 
culated scattering  with  Roberts  potential,  Or.    The  impact 
parameter,  b,  is  in  units  of  aQ  =  0.53(10"-°)  cm.    Columns  3-6 
represent  the  ratios  of  the  scattering  angles  calculated  with 
the  Bohr  potential  (B),  the  Csavinszky  potential  (C).  the 
I'iOliere  potential  (M),  and  the  Wedepohl  potential  (tf)  to  the 
corresponding  value  of  Or.    The  number  in  parentheses  after  an 
entry  is  the  power  of  10  to  be  multiplied  by  it. 


'■  "  ■  — 

b 

9b/qr 

9c/9R 

0m/9r 

»4-;/9r 

2.5000 

0.9937 fl) 

0.0558 

0.8581* 

1.0862 

0.90I49 

2-.U000 

0.1122 

0.0663 

O.8768 

1.0815 

0.8938 

2.3000 

0:1267 

0.0785 

O.S977 

1.0793 

0.3867 

2 ;2000 

0.1U38 

0.0933 

0.9170 

1.0725 

O.8767 

2.1000 

0.1635 

0.1113 

0.93U5 

1.0660 

0.3699 

2.0000 

0.1869 

0.1308 

0.91*89 

1.0571 

0.8618 

l.?000 

0.21U3 

0.1529 

0.9626 

1.01*81* 

0.851*6 

1.8000 

O.2I468 

0.1803 

0.971*5 

1.0368 

O.8I485 
0.8)430 

1.7000 

0.2856 

0.2100 

0.9837 

1.0253 

1.6000 

0.3321 

0.2UU6 

0.9903 

1.0129 

0.8333 

1.5000 

0:3385 

0.28U5 

0.991*3 

0.9991 

0.3350 

I.I4OOO 

0.14571 

0.3292 

0.9955 

0.9871 

0.8330 

1:3000 

0.5b 2 3 

0.3300 

0.9930 

0.9737 

0.8310 

1.2000 

0.61*81* 

0:1*360 

0.987U 

0.9622 

0.8305 

1:1000 

0.7823 

0.h977 

0.9789 

0:9527 

0.8313 

1.0000 

0.951*1* 

0-561*7 

O.967U 

0.91*50 

0.83.U5 

0.9000 

0.1179(1) 

0  .  6372 

0.95^5 

0.9)405 

0.8337 

0.3000 

0.11*81(1) 

0.7131 

0. 937U 

0.9397 

0.81*53 

0.7000 

0.189UU) 

0.7911 

0.9215 

0.9U30 

0.35)47 

0.6000 

0.2H82(1) 

0.8690 

O.9069 

0.9506 

0.6677 

0.5000 

0.3361(1) 

0-91*21 

0.8970 

0.9612 

0.8838 

0.1*750 

0.3614,8(1) 

0.9595 

0.3961 

0.96bU 

0.8885 

0.U500 

0.3971(1) 

0.9759 

0.895U 

0.9677 

0.3936 

0.1*250 

0.14339(1) 

0:9916 

0.3956 

C9710 

0.8989 

0. 14  000 

0.1*757(1) 

I.OO67 

0.896I* 

0.971*2 

0.90U5 

0.3750 

0.5238(1) 

1.0207 

0.3980 

0.9775 

o.<;io3 

0.3500 

0.579)4(1) 

1.0337 

0.9007 

0.9805 

0.9165 

0.3250 

0.61*U5(D 

1.01*51 

O.90I42 

0:9831* 

0.9229 

0 . 3000 

0.7211(1) 

1.055b 

0.9089 

0.986U 

0:9296 

0:2750 

0:8127(1) 

1.061*2 

O.9D46 

0:9886 

0.9363 

0.2500 

0.9235(1) 

1.0712 

O.9213 

0.9909 

0:91*32 

0.2250 

0:1060(2) 

1.0766 

0.9292 

0.9928 

0.9501 

0.2000 

0.1232(2) 

1.0802 

0.9333 

0. 991*1* 

0.9563 

0.1750 

0.11453(2) 

1.0812 

0.91*79 

0.9958 

0.9626 

TABLE  VIII  (continued) 


9B/»R 


9c/3R 


9M/Qr 


9,-j/QR 


0 . 1500 

0.l7li8(2) 

1.0802 

0.9583 

n  i  ?<o 

0.21^9(2) 

I.076I4 

0.9689 

0  1  000 

0.?76u( 2) 

1.0702 

0.979U 

n  ■  07  C  n 
u  •  U ( pu 

n  "37^8(2) 

0.5510(2) 

I.O607 

0.9888 

n  ocloo 

1.0U76 

0.9960 

0.C71),  f  2  ) 

I.OU63 

0.9965 

o  .0)1 

0.^876(2) 

1.0)452 

0.9968 

0.60)1^(2) 

1.0UU3 

0.9971 

0.6223(2) 

1-.OU32 

0.9973 

w  •  Wi4  c  \j 

0'.6)ilO(2) 

1.0U22 

0.9976  • 

o:o)iOc; 

\J  •  WIS  >JJ 

0.6606(2) 

1.0hl2 

0.9980 

n  rnon 

U  •  \Jjy\J 

0. 6813  ( 2) 

1.0U00 

0.9932 

0'  0T7£ 

0.70^0(2) 

1.0389 

0.998)4 

n  rnAo 

0:T?^9( 2) 
w  •  (  <-j  7  \  j 

1.0377 

0.9987 

0.7^00( 2) 

I.O365 

0.9989 

n  mm 

0'.77^  (2) 

1.0353 

0.9991 

ri-0^1  K 

\J  •  w  ?  J-  _5 

n.  8023(2) 

1.0 3U0 

0.9992 

0.8307(2) 

1 .0327 

0.99% 

0.8606(2) 

1.031U 

0.9996 

0'  0?70 

O.R922( 2) 

1.0299 

0.9997 

0.92^6(2) 

1.0285 

0-.9998 

U  •  Uc  14.U 

o.qf)C)P>(  2) 

1.0270 

1.0000 

0.0225 

0.9980(2) 
0.1037(3) 

1.0255 

1.0001 

0.0210 

r.0239 

1.0001 

0.0195 

0-.10?9(3) 

1.0222 

1.0002 

0-.01B0 

0;1122(3) 

1.0206 

1.0002 

0.0165 

0.1168(3) 

1.0189 

1.0002 

0.0150 

0.1216(3) 

1.0171 

1.0003 

0.0135 

0.1266(^) 

1.015U 

1.0003 

0.0120 

0.1319(3) 

1.0136 

1.0003 

0.0105 

0.1373(3) 

1.0118 

1.0002 

0.0090 

0.1U30(3) 

1.0100 

1.0002 

0.0075 

O.Ui88(3) 

1.0082 

1.0002 

0.0060 

0.15U9(3) 

1-.006U 

1.0002 

o.ooh5 

0.1610(3) 

1.00U7 

1.0001 

0.0030 

0'.  I67  3(3) 

1.0031 

.1.0001 

0.0015 

0.17 36(3) 

1.0015 

1.0000 

0.9969 

0.9979 
0.9991 
l.OOGU 
1.0016 
1.0017 
1.0017 
1.0018 
1.0018 
1.0019 
1.0020 
1.0020 
1.0020 
1-.0020 
1.0020 
1.0020 
1.0020 
1.0020 
1-.0020 
1-.0020 
1.0019 
1.0019 
1.0019 
1.0018 
1.0017 
1.0016 

1.0015 
1.001U 
1.0013 

1.0012 
1.0010 
1.0009 
1.0007 
1.0006 
1.000U 
1.0003 
1.0001 


0.9677 
0.9706 
0.9109 

0.9665 
0.9571 

0.9562 
0.9556 
0.95U9 
0.95U2 

0.9537 
0.9532 

0-.9527 
0.9522 
0.9518 
0.9515 
0.9512 
0.95H 
0.9512 
0.9513 
0.9516 
0.9521 
0.9529 
0.9538 
0.9550 
0.9565 
0.9583 
O.960I4 
0.9629 
0.9657 
0.9638 
0.9722 
0.9758 
0.9797 
0.9838 
0.9879 
0.9921 
0.9961 


TABLE  DC 


Calculated  differential  cross  sections  for  the  Bohr  potential  (B),  the 
Csavinszky  potential  (C),  the  Roberts  potential  (R),  the  Moliere  potential 
(M),  and  the  V/edepohl  potential  (W)  in  units  of       (aQ  =  0.53(10~°)cm) . 
The  values  of  impact  parameter,  b,  are  in  units  of  aQ.    E  =  h2.9  kev 
^Lab  =  100  kev),  Zj,  =  8,  and  Z?  =  6.    The  number  in  parentheses  after  an 
entry  is  the  pcwer  of  10  to  be  multiplied  by  it. 


2.500 

0.^615(9) 

0.^088(7) 

0.2592(7) 

2-.h00 

0.1711(9) 

0.2227(7) 

0.2211(7) 

2.300 

0.1076(9) 

0.1598(7) 

0.1598(7) 

2.200 

0.55^2(8) 

0.1119(7) 

0.1088(7) 

2.100 

0.2833(8) 

0.7932(6) 

0.7353(6) 

2.000 

0.1169(8) 

0-5532(6) 

0-5581(6) 

1.900 

0-75h7(7) 

0.38)46(6) 

0.3927(6) 

1.800 

0-1075(7) 

0.2623(6) 

0.2632(6) 

1.700 

0.20U6(7) 

0.1802(6) 

0.1826(6) 

1.600 

0.1033(7) 

0.12H(6) 

0.1229(6) 

1.500 

0:5225(6) 

0.8099(5) 

0.8137(5) 

l.hOO 

0.2577(6) 

0.5317(5) 

0.5259(5) 

1-.300 

0.1279(6) 

0.3U17(5) 

0.3301(5) 

1.200 

0.618U(5) 

0.2172(5) 

0.201)4(5) 

1.100 

0.2932(5) 
0.1362(5) 

0.1333(5) 

0.1211(5) 

1.000 

0.7973(h) 

0.6955(h) 

0.900 

0.6170(h) 

o.U558(b) 

0.3831(h) 
0.1998(h) 

0.800 

0.2695(b) 

0.2)467(14) 

0.700 

O.1127 (U) 

0.1237 (h) 

0.9795(3) 

0.600 

O.UU97(?) 

0.5636(3) 

0.hh66(3) 

o-.soo 

0.162H(3) 

0.2299(3) 

0.178)i(3) 

0.1475 

0.1201(3) 

0.1750(3) 

0.1352(3) 

0.U50 

0.9110(2) 

0.111h(3) 

0.10h8(3) 

0-.L25 

0.68U3(2) 

0.982)4(2) 

0.7959(2) 

o.koo 

0.50)48(2) 

0.7l2h(2) 

0.59?9(2) 

0.375 

0.1678(2) 

0.5118(2) 

0.h:,32(2) 

0.350 

0.26lil(2) 

0.1785(2) 

0.3179(2) 

0.325 

0.1863(2) 

0.261t0(2) 

0.2250(2) 

0.300 

0.1286(2) 

0.1797(2) 

0.1561(2) 

0.275 

0.8668(1) 

0.1190(2) 

0.1053(2) 

0.250 

0.5680(1) 

0.7593(1) 

0.6359(1) 

0.225 

0.3565(1) 

O.U67h(D 

0.h3H(l) 

0.200 

0.2165(1) 

0.27 37(1) 

0.2579(1) 

0.175 

0.1231(1) 

0.1508(1) 

o.lh52(l) 

0.150 

0.6508 

O.77H 

0.75h7 

0.125 

0.1101 

0.35U8 

0.3531 

0.100 

0.1278 

O.D4.15 

0.1U27 

0.2372(7 

0.l7  3h(7 
0.1357(7 
0.1003(7 
0.7319(6 
0.5279(6 
0.3800(6 
0.2719(6 
0.1897(6 
0.1318(6 
0.8889(5 
0.5832(5 

0.37hh(5 
0.23h3(5 
0.1h03(5 

0.30hh(h 

0.h335(h 

O.2256O4 

0.1076(U 

0.h7h7(3 

0: 1868 (3 

0-.lhl9(3 

0.1078(1 

0.8093(2 

0.6033(2 

0.hh29(2 

0.1201(2 

0.2267(2 

0.1565(2 

0.1055(2 

0.6881(1 

0.h321(l 

0.258h(l 

0.1h52(l) 

0.75hh 

0.3521 

o.lhlS 


0.3788(7 

0.2725(7 

0.2068(7 

0.151A  (7 

0.1132(7 

0.7990(6 

0.5572(6 

0.3290(6 

0.2710(6 

0.1822(6 

0.1190(6 

0.7710(5 

0.h33h(5 

0.2939(5 

0.1735(5 

0.9801(h 

0.5299(h 

0.2699(U 

0.1276(h 

0.5587(3 

0:2160(3 

0.1621(3 

0.1233(3 

0.9181(2 

0.6735(2 

0.h962(2 

0.35UK2 

0.2h79(2 

0.1692(2 

0.1127(2 

0.7300(1 

0-h5h2(l 

0.2697(1 

0.1511(1 

0.7398 

0.37U2 

0.15U9 


3ii 


TABLE  IX  (continued) 

b          4j(9)             6cO)              ^(e)  £M(S)  6w^8) 

0.075       O.L22?(-l)       o.U5U9(-D       0.U612(-1)  o.U578(-i)  o.$2U6(-l) 

0.050        0.9h^6(-2)        0.9887(-2)         0.1003(-D  0.9927(-2)  0.1232(-1) 
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TABLE  X 


Calculated  differential  cross  sections  for  the  Bohr*  potential  (B),  the 
Csavinszky  potential  (C),  the  Roberts  potential  (R).  the  Mbliere  potential 
(M),  and  the  Wedepohl  potential  (W)  in  units  of  a|  (aQ  =  0.53(10-°)  cm). 
The  values  of  impact  parameter,  b,  are  in  units  of  a  .    E  ■  21.1;  kev 
(Ej-v,  =  50  kev),  Zj_  ■  8,  and  Z2  =  6.    Tne  number  in  parentheses  after  an 
entry  is  the  power  of  10  to  be  multiplied  by  it. 


1 

b 

^  (0) 

A  (9) 

A  (9) 

A  (9) 

£  (9) 



2:500 

0.8773(8) 

0.7970(6) 

0.6713(6) 

0.6111(6) 

0.9597  (6 

2.!j.  00 

0.1,909(8) 

0.5716(6) 

0.5^17(6) 

O..Uh67  (6) 

0.68)42(6 

2.300 

0.2585(8) 

0.3970(6) 

0.^785(6) 

0.336U(6) 

0.5232(6 

2 . 200 

0.1317 ( 8) 

0.2812(6) 

0.27hO(6) 

0.250h(6) 

0.3851(6 

2 . 100 

0.1996(6) 

0.1968(6) 

0.l83h(6) 

0.2782(6 

2.000 

0.3778(7 ) 

0.139)j(6) 

0.1392(6) 

0.1326(6) 

0.2009(6 

I.900 

0.1926(7 ) 

0.9623(5) 

0.9786(5) 

0.9578(5) 

0.11418(6 

1.800 

0.9903(6) 

0. 661^(5) 

0.6757 (5) 

0.6816(5) 

0.9851(5 

1.700 

0.5166(6) 

0.lj6lU(5) 

0. Ii759(5) 

0.6776(5 

1  .600 

0.2598(6) 

0.30).°(5) 

n.  ?03K  ( £  1 

0.309h(5) 

0. 3300(5) 

0.1i51i)i(5 

1  con 

O'.l  *3in(M 

0.20)15(5) 

0.2227 (< ) 

O.oqRq^ 

0".6^in  (5) 

0.1338(5) 

0.1320(5) 

0".l!i70(5) 

0.1 932(5 

1.300 

fy.  3192(5) 

0.8635(h ) 

0.8312(h) 

0.9)i73(  1 ) 

O'.l  21h(5 

1  ■  ?oo 

0.1551  (5) 

0.512).  (). ) 

0.588)i()i) 

0'.7).29()i 

1.100 

0.739),()i) 

v>  •    J  J  ^  7  V  £4-  / 

0. 3058()i ) 

0.^^?7  ()i ) 
v  •       (  \  14  / 

0;)i?70()( 

X  .  \J\J\J 

O  -  POl  1  f)t  ^ 

0.1152(14) 

O  176T  ()i ) 

0*.?O^R( li 

O  ?)i770i 

0.900 

0.1559(h) 

0.970h(3) 

0'.1112(li) 

0.131)1(14 

0.800 

0.68)46(3) 

0.6238(3) 

0.5077(3) 

0.5720(3) 

0.6837(3 

0.700 

0.2877(3) 

0.31^9(3) 

0.21.96(3) 

0-.27UU(3) 

0.3252(3 

0.600 

0.1152(3) 

0.1U63(3) 

0.11U0(3) 

0.1219(3) 

0.  II4  30(3 

0.500 

0.U222(2) 

0.5867(2) 

O.U601(2) 

O.U803(2) 

0.5556(2 

0.1*75 

0.3152(2) 

0-.UUUK2) 

0.3508(2) 

0.3628(2) 

0.14166(2 

0.I450 

0.2386(2) 

0.338h(2) 

0.2703(2) 

0.2777(2) 

0.3171(2 

0.I425 

0.1783(2) 

0.25U6(2) 

0.2056(2) 

0;2102(2) 

0.238)4(2 

o.uoo 

0.1320(2) 

0.189M2) 

0;l5h3(2) 

0.1569(2) 

0.1767(2 

0.375 

0.9653(1) 

0.1381(2) 

0.11kl(2) 

0.1156(2) 

0.1290(2 

0.350 

0.6955(1) 

0.9876(1) 

0.8279(1) 

0.8371(1) 

0.92142(1 

0-.325 

0.1.927(1) 

0.6920(1) 

0.5895(1) 

0.59^0(1) 

0.6149^(1 

0.300 

0.3U16(1) 

0;h7°0(l) 

0.U102(1) 

0.hl22(l) 

0;U)i68(l 

0.275 

0.2315(1) 

0.3lii6(l) 

0.2778(1) 

0.2790(1) 

0.2991(1 

0.250 

0.152^»(1) 

0.2025(1) 

0.l82h(l) 

0.1828(1) 

0.19h0(l 

0.225 

O.9667 

0.125U(1) 

0.1153(1) 

0.1155(1) 

0.1215(1 

0.200 

0.5885 

0.7U19 

0.6950 

O.6962 

O.728I4 

0.175 

0.3387 

0.U1L3 

0.3951* 

0.3959 

O.I4I26 

0.150 

0.1813 

0.211*9 

0.2087 

0.2089 

0.2179 

0.125 

0.6829(-L) 

0.1013 

0.9988M) 

0.9979^-D 

0.10U9 
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TABLE  X  (continued) 


b 

46(o) 

6Q(Q) 

^(9) 

0.100 

0.075 
0.050 

0.3768( 
0.13li8( 
O.3601( 

-1) 
-1) 
-2) 

O.U187(-D 
o.iU5o(-l) 
0.^7U3(-2) 

0.hl78(-l) 
0.1h57  (-D 
0.3769(-2) 

o.)4l65(-D 
o.il,5K-i) 
0.3750(-2) 

0.1i)jh9(-D 
0.1585(-D 
0.]4l68(-2) 
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where  f(t^"'^)  is  a  scaling  factor  called  the  reduced  cross  section.  The 
factor  t"3/2  was  introduced  so  that  for  a  potential  proportional  to  r~^, 
f(t1/2) 

is  constant.    In  the  present  calculations, 

6  (9)  =  $f  =  —-^ts  (27) 
is  evaluated  by  equations  (2li).    From  equation  (2j>) 

Substituting  the  above  equation  into  equation  (27)  we  get 

d/=  hn  £(Q) 
dt  £2 

and  f(t^/^)  is  given  by  the  relation 

/ 

f(t1/2)  =  8tV2_£i9l  (28) 

Except  for  potentials  of  the  form,  V(r)  =  c/rn,  f(t1/2)  must  be  evaluated 
numerically.    The  calculated  values  of  0  and  £(Q)  were  substituted  into 
equations  (25)  and  (27)  to  evaluate  f(t1/2).    Results  of  these  calculations 
are  presented  in  figures  1  and  2  and  in  tables  XI  through  XV.    These  curves 
can  be  used  for  any  combination  of  Z-^,        and  E.* 

This  completes  a  detailed  comparison  of  scattering  angle  and  differential 
cross  section  for  various  forms  of  the  interatomic  potentials.    The  discus- 
sion and  conclusion  are  in  the  next  section. 


"VJhere  the  classical  calculations  are  valid. 
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TABLE  XI 


Calculated  values  of  tV?  and  reduced  cross  section,  f (tV^) t 
for  the  Bohr  potential,    t  and  f(tV')  are  defined  in 
equations  (25)  and  (28).    The  number  in  parentheses  after  an 
entry  is  the  paver  of  10  to  be  multiplied  by  it. 


tl/2  f(tl/2)  tl/2  f(tl/2) 


0.  £o5>05  V 

-3; 

U  .J.i  i3oO\ 

-2 ) 

U  .11  v?m  i; 

0«  5$.5oo  v 

-3; 

-2J 

U m 5  (  0 

O.U7-T/5  I 

-3  J 

U .  ( £0)y { 

-d) 

,0  1  G'7 1 ,  £  (  1  \ 

0  .C£  -<lo  v 

-3; 

-2,) 

U • cy^oc \1 ) 

O.oUoo^ ( 

-3) 

(J.lJ.O  _V(, 

-I; 

r\       ry->  r1  L  f  ~\  \ 
V.  _j\J(  'pOK  L) 

O.lljoOv, 

-2) 

0 .130511 
o.ioiyii 

-l; 

0.  .169311; 

0.13579 

0.1552U( 

-2) 

-i) 

0.32003V1J 

O.I3-I-7O 

0.?0oo5( 

-2) 

0.2  jo50i 

-1) 

0  ♦  337  32  (.1; 

O.12017 

0 . 2cG65 ( 

-2) 

0.291U3', 

-t  \ 

-l; 

0ou^u3l.l ) 

0.12[ju2 

0.3c2!i2(, 

-2) 

0.37U15( 

-1) 

C3&02U(j.) 

0.1207U 

i  <-?  0  Q  f 

O.5I798I 

-2) 

O.U6777* 

-1) 

0.37278(1; 

0.11730 

O.70oo9v 

-2  J 

0.59i3Li  1 

-  N. 
-1J 

0.38610(1) 

0 .11338 

0. 96875  ( 

-2) 

0.7U89UI 

-1) 

0..u00  35(lj 

0.10991 

0.1328U( 

-1) 

0 . 93  37  9 ( 

-1) 

0.)4l5h8(l) 

0.106U0 

0.18 313 ( 

-1) 

0.11600 

0.1i3169(l) 

0.102146 

0. 25389 ( 

-1) 

0.1ui6U 

o.)iU905(D 

0.987  97  ( 

-1) 

u./55kol 

-1J 

U .1/ Ooy 

'J.UD^OOilJ 

0.951o7( 

-i; 

O.U9825( 

-1) 

0.21U71 

0.[|R?63(D 

0.91660( 

-1) 

0.7088U( 

-1) 

0.25650 

0.5090U(1) 

O.R8070( 

-1) 

0.10225 

0.30972 

0.53209(1) 

0.8l4237( 

-1) 

O.15063 

0.35U36 

0.55697(1) 

0.808hO( 

-1) 

0.16659 

0.35775 

0.58363(1) 

0.773hU( 

-1) 

0.18h76 

0.37093 

•  0.612U2(1) 

0.73796( 

-1) 

0.20528 

0.38135 

0.6U^31(1) 

0.70509  ( 

-1) 

0.22870 

0.38897 

O.676UMI) 

O.67 095 ( 

-l) 

O.2556O 

0.39569 

0.71192(1) 

0.63368( 

-1) 

0.28666 

0.I10070 

0.71x958(1) 

0. 60866 ( 

-1) 

0.32283 

O.U0387 

0.78925(1) 

0.57680( 

-1) 

0.^6529 

0.^0386 

0.83063(1) 

0.5505K 

-1) 

O.U1582 

O.UOII18 

0.87^05(1) 

0.5'2l|28( 

-1) 

0.U765U 

0.39597 

0.91557(1) 
0.95676(1) 

0.50068( 

-1) 

0.55077 

0.38583 

O.U797l( 

-1) 

0.6U338 

0.371U5 

0.99U92(1) 

O.U6l7U( 

-1) 

O.76I76 

0.35115 

0-10279(2) 

O.Uj72U( 

-1) 

0.9185U 

0.32U92 

0.10536(2) 
0.10700(2) 

O.U358U( 
O.U2953( 

-1) 
-1) 

TABLE  XII 


Calculated  values  of  and  reduced  cross  section.  f(tV2), 

for  the  Csavinszky  potential,    t  and  f(tV2)  are  defined 
in  equations  (25)  and  (28).    The  number  in  parentheses 
after  an  entry  is  the  power  of  10  to  be  multiplied  by  it. 


t1/2 

f(tV2) 

t1/2 

f(tV2) 

o.uooiM 

-2) 

0.127UU 

0.10231(1) 

0.2hli82 

O.U6285I 

-2) 

0.1L(225 

0.13361(1) 

O.217I4O 

0.535'20( 

-2) 

0-1578U 

0.181*96(1) 

0.1851*3 

0.61910( 

-2) 

0.17U10 

0.28291(1) 

0.11*1*23 

O.7I6901 

-2) 

0.18828 

0.291*69(1) 

0.12909 

0-.8335U' 

-2) 

0.20638 

0.30)4  00(1) 

0.12)453 

0. 96898 ( 

-2) 

0.22$kh 

0.31390(1) 

0.12173 

0.11302( 

-1) 

0.2LUU2 

0-321*38(1) 

0.11810 

0.132101 

-1) 

0.26762  ( 

0.33556(1) 

0.115)42 

0.l5U6l( 

-1) 

0.26881 

0.31*731(1) 

0.1121*5 

0.18173( 

-1) 

0.3131U 

0.35983(1) 

0.10911 

0.21L|00( 

-1) 

0.33569 

0.37320(1) 

0.10605 

0.25^7  ( 

-1) 

0.35851 

0.387lib(l) 

0.1027 1* 

0.30lLi7( 

-1) 

0. 38310* 

0.1*0262(1) 

0. 996)4)4  < 

-1) 

0O6C58I 

-1) 

o.ho)4ou 

0.1*1386(1) 

0.96170( 

-1) 

0-U3U95( 

-1) 

042265 

0.1*3630(1) 

0.93023( 

-1) 

0.52982( 

-1) 

0.U367U 

0.1*51*92(1) 

0.89930( 

-1) 

0.65UU9C 

-1) 

0.W455  2 

0.1*71*95(1) 

0. 86)481  ( 

-1) 

0. 82353 ( 

-1) 

0.liU522 

0.1*9651(1) 

0.8?205( 

-1) 

0.106)48 

0  .Ltlj.221 

0.51971(1) 
o.5UU69(D 

0. 80011 1 

-1) 

0.1U279 

0.U3115 

0-767)45( 

-1) 

0.15H95 

0.U1932 

0.57162(1) 

0.73533( 

-1) 

O.I6863 

o.U059^ 

0.60062(1) 

0.70)»19( 

-1) 

0-.18UU6 

0.39725 

0.63181(1) 

0.67197 ( 

-1) 

0'.  20266 

0.39270 

0.66539(1) 

0.6hl29( 

-1) 

0.22368 

0-.383U1 

0.70128(1) 

0.61121( 

-1) 

0.2U852 

0.37U30 

0.73958(1) 

0.58172( 

-1) 

0-.27781 

0.36)479 

0.7800)4(1) 

0. 55 3 86 ( 

-1) 

0.31307 

0.35529 

0.822?3(1) 

0.52697 ( 

-1) 

0.35561 

0.31|U70 

0.8658U(1) 

0.50223( 

-1) 

o.)j08ou 

0.3^232 

0.90951(1) 

O.U7978( 

-1) 

0.1*7363 

0.31996 

0.95205(1) 

0.l4$9U2( 

-1) 

0.55716 

0.30LI93 

0.99157(1) 

O.Ulj?06( 

-1) 

0.66673 

0.2878U 

0.10259(2) 

0.1*2807  ( 

-1) 

O.8U469 

0.26860 

0.10526(2) 
0.10697(2) 

0.1*17  01  ( 
0.hl067( 

-1) 
-1) 

TABLE  XIII 


Calculated  values  of  and  reduced  cross  section, 

f(t1/2)>  for  the  Roberts  potential,    t  and  fit1/2)  are 
defined  in  equations  (?5)  and  (28).    The  number  in 
parentheses  after  an  entry  is  the  power  of  10  to  be 
multiplied  by  it. 


tV2 


f(tV2) 


tV2 


f(tV2) 


0-U6659 
0.52718 
0.59UU1 
0.67587 
0.76807 
0.878*1 
0.10055 
0.11589 
0.13h2^ 
0.15625 
0.18278 
0.21512 
0.25510 
0.3052)4 
0.36826 
O.iiJ-968 
Oo56l7 
0.6938U 

0.89537  (-D 

0.117U9 

0.159U3 

0.17  312 

O.I8863 

0.2061? 

0.22626 

0.2U939 

0.27613 

0.307U3 

0.38873 
0.bh260 
0.50910 
0:59239 
0.70153 
0.8h700 


(-2) 
(-2) 
(-2) 
(-2) 
(-2) 
(-2) 
(-D 
(-1) 
(-1) 
(-D 
(-D 
(-1) 
(-D 
(-1) 
(-D 
(-1) 
(-D 
(-D 


0.16961 
0.21077 
O.20U61 

0.216U5 

0.22922 

0.2U?6l 

0.25712 

0.26387 

0.28!45b; 

0.30207 

0.32009 

0.^3725 

0.3532)4 

O.37272 

0.39038 

0.ij07h5 
0.U2158 

O.U3935 
0.U5293 
O46658 
0.b6l26 
O.U5188 
0.h529U 
O.UU936 
0.!^320 

0-.-U379U 
O.U3H3 
0.U211U 
O.I4II27 
0. 3981*2 
0.33313 
0.^661i7 
0.?h623 
0.^05 

0.295U5 


0.10516(1 
0.13571(1 
0.13597(1 

0.282)47(1 
0.29)407(1 
0.30328(1 
0.^1310(1 
0.32^6(1 
0.33]4).9(1 
0.3U6l7(l 
0.^5856(1 
0.37180(1 

0.38597(1 
0.)4  0102(1 
O.U72l(l 
O.U3)jU6(1 

o.U5 305(1 

O.U7297(l 
0.U9UUK1 
0-5l75U(l 
o.5U2UU(i 
0.5693UU 
0.5983)4(1 
0.62956(1 
0.66311(1 
0.69911(1 
0-737U3(l 
0.7780)4(1 
0.820)49(1 
0.86)tl6(l 
0.90813(1 
0.95097(1 
0.99078(1 
0.1025h(2 
0.l052)j(2 
0.10697(2 


0.26)450 

0.22977 

0.19111 

O.H4561 

0.1*071 

0.12)470 

0.12218 

0 . 11905 

O.II572 

0.11283 

0.10955 

O.IO576 

0.10269 

0-9935 2 ( 

-1) 

0. 96115 1 

-1) 

0.928)4  2  ( 

-1) 

0.89U77( 

-1) 

0. 86301 ( 

-1) 

0.82835( 

-1) 

0.79577 ( 

-1) 

O.76297  ( 

-1) 

0.72983C 

-1) 

0.69830( 

-1) 

0.667U2( 

-1) 

0.63613( 

-1) 

0.606h7( 

-1) 

0.57706( 

-1) 

0.5)4885  ( 

-1) 

0.5228h( 

-1) 

0- It  97  35  ( 

-1) 

0-.h7lj89( 

-1) 

0.lj$li93( 

-1) 

0.)43779( 

-1) 

0.h2373( 

-1) 

0.)4128h( 

-1) 

0.U0662( 

-1) 

TABLE  XIV 


Calculated  values  of  tV2  and  reduced  cross  section, 
f(t-L/^),  for  the  >ioliere  potential,    t  and  f(tV2) 

are 

defined  in  equations  (25)  and  (28).    The  number  in 
parentheses  after  an  entry  is  the  power  of  10  to  be 
multiplied  by  it. 


tV2 

f(tV2) 

0'.5062O 

:-2) 

0.1982U 

0.570h7( 

-2) 

0.20738 

0.6h298( 

-2) 

0.23236 

0-.7  23951 

-2) 

0.2h5l8 

0.81906( 

-2) 

0.25908 

0.92757< 

-2) 

O.27II4O 

0.105 61 ( 

-1) 

0.288^6 

0. 12027 ( 

-1) 

0.^0)467, 

0.13767( 

-1) 

0.31883 

0.15812( 

-1) 

0.-(356Li 

0.182U6( 

-1) 

0.3U782 

0.21210( 

-1) 

0.35850 

0.2h85K 

-1) 

0  .  3701U 

0.29 376 ( 

-1) 

0.38266 

0.3507 3 ( 

-1) 

0.39121 

O.U2)i75( 

-1) 

0.39710 

0.^08( 

-1) 

0.U0h30 

0. 6566)4  ( 
0.8U396( 

-1) 

0.U11U7 

-1) 

0.U1660 

0.1117U 

O.U2658 

0.15^28 

O.U3337 

0.16? 07 

0-U2613 

0.18256 

O.U2255 

0.20028 

O.U1912 

0.22059 

O.U1718 

0.2U39li 

O.UD4I9 

0.27102 

O.U1082 

0.30259 

O.UOh63 

o.3hooi 

0.39633 

0.38h72 

0.38698 

O.h3o90 

0.37b 90 

0.50571 

0.3600b 

0.53983 

0.3U162 

O.69877 

0.^1925 

0.8UU89 

0.29n2 

4.1/? 

f(tV2) 

0.10500(1) 

0.26255 

0.13573(1) 

0.228U6 

0.19070 

0.2o3U5 (l) 

O.D4563 

0.29509(1) 

0.12685 

0.^0hii5U) 

0.12688 

r\    nl  or*  f  ~\  \ 

O.^lij  25(1) 

0.12298 

0.  ^2)|Dy(i; 

0.11887 

0.11586 

0-3U755(l) 

0.11317 

0 . <.o002(l ) 

0.10958 

0.37335(1) 

0.106U7 

■ 

0.387U9(1) 

0.10310 

0.140269(1) 

0.9967K 

-1) 

O.U18o9(l) 

0. 96611 ( 

-1) 

O.U362Lt  (1) 

0.93253( 

-1) 

0.)45U86(1) 

0. 89817 ( 

-1) 

0.)j7U90(1) 

0. 86606 ( 

-1) 

0.]i9oU0(l) 

0.33308( 

-1) 

0.51959(1) 

0.79979( 

-1) 

o.5Wi55(D 

0.76836( 

-1) 

0.5711414(1) 

0.73U56( 

-1) 

0.60050(1) 

0-70273( 

-1) 

0.63170(1) 

O.67 250 ( 

-1) 

0.66523(1) 

0.6b06b( 

-1) 

0.70121(1) 

0.61037 ( 

-1) 

0.73950(1) 

0.5 8160 ( 

-1) 

o.77997(D 

o.5573U( 

-1) 

0.82230(1) 

0.5269^ 

-1) 

0.86578(1) 

0.50226( 

-1) 

0.90950(1) 

O.L79b7( 

-1) 

0.9520?(1) 

o.)45929( 

-1) 

0.99155(D 

0.bbl93( 

-1) 

0.10259(2) 

0.U2790( 

-1) 

0.10526(2) 

0.hl708( 

-1) 

0.10697(2) 

0.14107  9  ( 

-1) 

TABLE  XV 


Calculated  values  of  tV2  and  reduced  cross  section, 
f(tV2),  for  the  VJedepohl  potential,    t  and  f(tV2)  are 
defined  in  equations  (25)  and  (28).    The  number  in 
parentheses  after  an  entry  is  the  power  of  13  to  be 
multiplied  by  it. 


tl/2  f(tl/2)  tl/2  f(ti/2) 


O.U2207 (-2) 
0-.U713K-2) 
0.528h9(-2) 
0.593U2(-2) 
0.66836(-2) 

0.755)41  (-2) 
0.86ll0(-2) 
0-9852M-2) 
0.1132U(-1) 
0. 13085 (-1) 
0.152U3(-1) 
0.17903(-1) 

0. 21175 (-1) 

0.25335(-D 

0.30616(-1) 

0.^7500(-l) 

0.U6633(-1) 

0.59033(-l) 

0.76U7UC-D 

0-.  1019b 

0.1h089 

0.15385 

0.168J2 

0.18533 

0.20)477 

0.2271U 

0.25319 

0.2839b 

0.320U6 

0.36U55 

0-U1810 

O.U8U33 

0.56799 

0.67627 

0. 82058 


0.183U6 

0.18381 

0.19662 

0.20748 

0:21773 

0.22188 

0.22919 

0.2U51 

0.2535U' 

0.26290 

0.27155 

0.28502 

0.29565 

0.30791 

0.32079 

0.33296 

0.3L621 

0.35765 

0.3676U 

0.38129 

0.38918 

0.33026 

0:38005 

0.376)48 

0-37531 

0.37u57 

0.37026 

0.36552 

0.35380 

O.35179 

0;  3)4^7  2 

0.332)16 

0.318J35 

0.30139 

0.23112 


0.10206(1) 
0.13138(1) 
0.17811(1) 
0.26358(1) 
0.27361(1) 
0.23156(1) 
0.2899U(1) 
0.29883(1) 
0.30820(1) 
0.31816(1) 
0.32866(1) 
0.33991(1) 
0.35-180(1) 
0.36)458(1) 
0.37813(1) 
0.3927HU) 
O.U08.U2(l) 
O.)42530(l) 
0.1^31*3(1) 
O.I4631KI) 

o.U8U56(l) 
0.50781(1) 
0.53310(1) 
0.56085(1) 
0.59119(1) 
0.621.50(1) 
0.66106(1) 
0.7011U(1) 
0.7UU92(1) 
0.79230(1) 
0.8U28K1) 
0.895)47(1) 
0.9)480)1(1) 
0.99720(1) 
0.10382(2) 
0.10658(2) 


0.25626 

0.22626 

0.19093 

O.D4527 

0.12953 

0.1267U 

0.122U7 

0.1190U 

0.115)42 

0.11218 

0.10837 

0.10U73 

0.10101 

o.97U96(-l) 

0.93839(-D 

0.897U5(-D 

0.85989(-D 

0.82U36(-1) 

0.78669(-l) 

0.7UU10(-1) 

0.7O55U(-D 
0.66970(-l) 

0. 62995  (-D 
0-.59118(-l) 
0.553U8(-1) 
o.5i552(-D 
O.U786U(-l) 
O.Uh2!i7(-l) 
0.1i0796(-l) 
0.3755U(-D 
0:31iU6l(-l) 

0.31670(-1) 
0.2923!t('-l) 
0.27202f-l) 
0. 25600 f-1) 
0.2U658(-1) 


V.    DISCUSSION  AND  CONCLUSION 

The  procedure  used  in  the  present  calculations  was  checked  by  a 

8  i 

comparison  of  the  present  calculations  to  those  of  Everhart  et  al.  We 
have  extended  the  above  mentioned  calculations  to  include  small  scattering 
angles. 

Values  of  £(9)  and  Q(b)  were  calculated  for  the  Bohr,  Moliere, 
Csavinszky,  Roberts,  and  Wedepohl  potentials.    Table  VIII  illustrates 
the  sensitivity  of  the  scattering  angle  for  a  wide  range  of  values  of 

0 

impact  parameter,  b.    We  note  that  for  b^l  A,  the  scattering  angle  from  an 
interaction  of  a  50  kev  oxygen  ion  with  a  carbon  target  atom  with  the  Bohr 
potential  is  about  20  times  smaller  than  angles  with  other  potentials  con- 

o 

sidered  in  this  work.  For  smaller  values  of  b  (^0..3  A)  all  the  potentials 
predict  the  same  scattering  angles  within  about  10  per  cent. 

Csavinszky  reports  that  the  potential,  which  he  describes  in  a  closed 
form  (equation  (19)),  is  in  reasonable  agreement  with  the  Thomas -Fermi 
potential  up  to  intermediate  distances.    However,  the  Csavinszky  potential 
has  been  claimed  to  have  more  desirable  behavior  at  large  distances.  The 
Firsov  form  of  the  interatomic  potential  was  used  in  all  calculations 
involving  the  Csavinszky,  Roberts,  and  Wedepohl  potentials.    The  difference 
in  the  screening  length,  as  defined  by  Thomas-Fermi,  and  that  of  Firsov 
differs  by  less  than  H  per  cent. 

The  reduced  cross  sections,  f(t^/^),  for  the  different  potentials  are 

given  in  figures  1  and  2  to  provide  a  comparison  with  the  results  of  the 

Thomas-Fermi  potential  obtained  by  Lindhard  et  al.  5    It  is  noted  that  the 

reduced  cross  section  for  the  Moliere  potential  agrees  reasonably  well  with 

1  /2 

the  Thomas-Fermi  values  over  a  wide  range  of  t  '  . 


U6 

The  classical  calculations  are  generally  valid  when  0  is  greater  than 
0*       X/2rca.    When  the  classical  results  may  be  in  error,  <5"(0)  can  be 
calculated  with  the  Born  approximation.    The  two  methods  complement  each 
other. 

The  present  calculations  can  be  very  conveniently  used  in  Monte  Carlo 
calculations  of  elastic  energy  loss  for  ions  in  amorphous  targets.  Some 
preliminary  work  has  already  been  completed  and  further  investigations  are 
planned. 
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ABSTRACT 

Several  screened  potentials  are  used  in  the  calculation  of  the 
differential  cross  section  and  the  scattering  angle  as  a  function  of  impact 
parameter.    The  numerical  results  of  Everhard  et  al.     (for  the  Bohr 
potential)  pre  extended  for  small  angles.    The  reduced  differential  cross 
section  is  presented  as  a  function  of  the  dimensionless  variable  t1/2  for 
five  potentials.    These  results  can  be  used  for  any  Zj_,  Zp,  and  energy 
where  the  quantum  effects  are  negligible.    The  calculations  were  performed 
on  the  IBM  360/50  computer  at  Kansas  State  University.    The  basic  computer 
program  can  be  used  in  the  Monte  Carlo  calculations  of  elastic  energy  loss 
for  ions  in  matter. 


